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■ ABSTRACT 

(N 

Context. We carry out a detailed 2D study of the ionised gas in the local universe galaxy UGC9837. In nearby galaxies, like the galaxy 
in question here, the spatial distribution of the physical properties can be studied in detail, providing benchmarks for galaxy formation 
models. 

r ^ ' Aims. Our aim is to derive detailed and spatially resolved physical properties of the ionised gas of UGC 9837. In addition, we derive 

I an integrated spectrum of the galaxy and study how varying spatial coverage affects the derived integrated properties. We also study 

<~i , how the same properties would be seen if the galaxy was placed at a higher redshift and observed as part of one of the high-z surveys. 

CL(' Methods. UGC9837 was observed using the PMAS PPAK integral field unit. The spectra are reduced and calibrated and the stellar 

and ionised components separated. Using strong emission line ratios of the ionised gas, the source of ionisation, the dust extinction, 
the star formation rate, the electron density and the oxygen abundance derived from a total integrated spectrum, central integrated 
spectrum and individual fibre spectra are studied. Finally, the same properties are studied in a spectrum whose spatial resolution is 
^ ■ degraded to simulate high-z observations. 

Results. The spatial distribution of the ionised gas properties is consistent with inside-out growing scenario of galaxies. We also 
find that lack of spatial coverage would bias the results derived from the integrated spectrum leading, e.g., to an under-estimation of 
ionisation and over-estimation of metallicity, if only the centre of the galaxy was covered by the spectrum. Our simulation of high-z 
^ ■ observations shows that part of the spatial information, such as dust and SFR distribution would be lost while shallower gradients in 

ly-^ ' metallicity and ionisation strength would be detected. 

Key words. Galaxies: individual: UGC9837 - Galaxies: abundances - Galaxies: ISM Galaxies: stellar content ~ Techniques: spec- 
troscopic - Stars: formation 
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1. Introduction galaxies (eg.. iRosales-Ortega et al]|2010t ISanchez et all 1201 lb . 

^, , , . , . ^ , , , . , Previous works have explored the use of different integral 

The local universe galaxies are fundamental anchor points when ^^j^ spectrophotometers for a detailed study of nearby galax- 

studying the evolution of galaxies with cosmologica time Their -^^ particular, the SAU RON project ( de Zeeuw et al.lM) 



apparent scale-lengths are small enough to derive relatively high .^^ extension Adas3D (ICappellari et all 1201 Ih are focused 

resolution information on their physical properties with the cur- ^^^^ early-type galaxies and bulges of spirals at 

H , rent technology. This information, in turn, provides constraints ^ ^ q ^^^^ ^^^-^ kinematics and stellar- populations. 

■ *at any theory of galaxy evolution need to meet. Th^ j^^^ ^j^^ ^^^^ ^^^.^ ^^.^^^^ ^ 42 Mpc) and 

Cold Dark Matter paradigm (ACDM) is currendy the most ^j^^ field-of-view (FOV) covered by the SAURON insti-ument 

widely accepted model describing the formation of structures (33,, ^ 4^,,^^ ^j^^^ ^^^^ ^^^^ restricted to the central part 

m the Universe In case of the spii-al galaxies, this niodel leads j^^.^^ ^^^^ ^^^^ j^^l j^^j^^ ^^^^ 

to t he so-called mside-out formation (see, ^^^^^^^ ^ ^j^^ ^^^q^ yiRUS-P spectrograph 

e.g.| White & Fren]^ [MIl ^^^ fcpr|1199a, i.e. the ^33 ^^J -^^^ j.^i,,^ B^l,„, 201O) , the DiskMass 

galaxy bulge forms first and is followed by the more slow pro- ^ combining PPAK and th e Sparse Pak spectrograph (146 



cess of disc formation. This scenario leads to various predictions ^^^^.j j^^.^^^ Bersha dT^M. and the PINGS 



about the galaxy structure such as negative age gradients, nega- ^^^^^^ ^^^^^^ ^^^j^^ galaxies, |Rosaies-Ortega eialj|20ld). 



tive gradients of the gas-phase heavy element abundances, and ^j^^ pp^^ particular, the mosaicking designed 
decreasing amount of dust with increasing galactocentric dis- pj^^^ ^ ^jj^^^^ j^^^ .^^^ 
tance. These predictions can be tested observationally ^j^^ ^j^^j^ extension of the galaxies and to explore the two- 
Integral field spectroscopy provides a powerfu tool for dimensional metallicity structure of disks. The galaxies of our 
studying the distribution of physical properties in well-resolved 
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sample (iMarmol-Oueralto et al.l201ll) were selected in size to fit 
in the FOV of the PPAK instrument, and therefore, to map their 
physical properties in their whole extension in just one pointing. 

We present here a detailed study of ionised gas in one galaxy, 
UGC9837. This Local Universe (z ~ 0.009) galaxy is posi- 
tioned nearly face-on towards us, making it a good candidate 
for 2D galactic disk studies. The data were observed as part of 
the PMAS/PPAK, one of the wo rlds widest FOV IFU, survey o f 
48 nearby galaxies presented bv lMarmol-Oueralto et alJ ( 1201 ll) . 
This su rvey was part of the preparatory phase of the CALIFA 
survey (ISanchez et al.| [2010; Sanchez et al. 2011), a much larger 
IFS survey comprising observations of 600 nearby galaxies. 

We use our data to derive for UGC9837 2-D distributions 
of: 1) local nebular reddening estimates based on the Balmer 
decrement; 2) oxygen abundance distributions based on a suite 
of strong line diagnostics incorporating reddening-corrected He, 
Hy6, [O ii], [O III], [N ii], and [S ii] line ratios; 3) measurements of 
ionisation structure in H ii regions and diffuse ionised gas using 
the well-known and most updated forbidden-line diagnostics in 
the oxygen and nitrogen lines; 4) the distribution of the surface 
star formation rate (SFR) traced by the Ha line intensity. The 
derived spatial distribution of these properties allows us to test 
the inside-out growing scenario of galaxies. 

In addition to the 2-D maps, we also study the integrated 
values of the above mentioned properties. The advantage of 
integral-field spectroscopy in this respect is that the observed 
spectra can be combined to produce an integrated spectrum of 
the object, using the IFU as a large aperture spectrograph. From 
this high signal-to-noise integrated spectrum the real average 
spectroscopic properties of the galaxy can be derived. This way 
no assumptions need to be made, in contrary to previous stud- 
ies (based on long-slit spectroscopy, for example) attempting 
to describe the average spectroscopic properties of the galaxy 
based on individual spectra taken at different regions. This way 
the integrated properties of the galaxies in the Local Universe 
can be compared with the corresponding properties of the unre- 
solved high-redshift objects. The most similar approach would 
be the spectrum derived by using drift-scanning technique (e.g. 
[Moustakas & Kennicutt 2006a). The advantage of using an IFU 
with respect to the drift-scan technique is that all the spectra are 
obtained simultaneously and that the IFU spectra allows a com- 
parison between the integrated and the spatially resolved prop- 
erties of the galaxy. 

Finally, based on our Local universe data, we study how 
the spatial properties of ionised gas could be seen in obser- 
vations at higher redshift, such as those carried out by SINS 
(Forster Schreiber et al. 2006) or MASSIV (Epinatetal. 200§ 
fOuevrel et al.ll2009 ) surveys with VLT/SINFONI. 

The article is structured in the following way: Sec. 2 presents 
the galaxy UGC9837; Sec 3 describes the observations and data 
reduction; in Sec. 4 the data analysis is presented; in Sec. 5 the 
results are discussed and in Sec. 6 the summary and conclusion 
are given. Throughout the paper we adopt a standard ACDM 
cosmology with Ho = 70.4, Q,„ = 0.268 and = 0.732. 

2. UGC9837 

UGC9837 is a spiral galaxy of magnitude 14.6 in the g band, 
located at r edshift z ~ 0.009, an d morphologically classified 
as SAB(s)c ( de Vaucouleurs et al .''T991). In addition to the op- 
tical (SDSS), it has pr eviously been observed at other wave- 
lengths: UV (GALEX , |Gilde Paz et al. 2007), NIR (2MASS, 
ISkrutskie etal] l2006h and (IRAS, .Beichman et alJ 1 19881) . see 




Fig. 1. SDSS (r band) image of UGC 9837. The cyan hexagon 
shows the approximate position of the PPAK pointing. The di- 
ameter of the green circle corresponds to the apparent diameter 
quoted in HyperLeda. North is up and East is left. 



Table [T] In additi on, this galaxy was analysed as part of the 
DiskMass Survey (iBershadv et al.ll2010h . 

T he distance to this galaxy has been estimated to be 41.2 
Mpc (lTullvlll988h with a distance modulus of 33.07 mag. The 
galaxy is extended (2.0' x 1.8', N ilsonlfT973 ^ and nearly face-on 
which permits the spatial properties of the disk to be well stud- 
ied. The total apparent corrected magnitude according to LEDA 
database is 1 3.65 mag, thus giving an absolute magnitude of 
- 19.42 mag. lPohlen & Truii llo (2006) reported the existence of 
a bar of length R <9" (~ 1 .64 kpc at the distance of the galaxy). 
These authors also found a break in the surface brightness profile 
at R ~ 50" , which almost corresponds to the optical diameter, 
D25 - 2'. 

3. Observations and data reduction 

UGC9837 was observed on 23rd June 2007. The night was not 
photometric and the extinction was high (Ay > 0.2) and vari- 
able. During the night the seeing varied from 1 to 1.8". The ob- 
servation were carried out using the 3.5 m telescope of the Calar 
Alto Observatory and the Potsdam Multi Aperture Spectr ograph, 
PMAS (fe oth et al. 2005) in the PPAK mode ( Verheiie netan 
!2004tlKelzet al.ll2006l) . 

Fig. 1 shows the SDSS r-band image of UGC 9837 with 
the approximate position of the PPAK FOV over-imposed. As it 
can be seen, a single PPAK exposure covers approximately the 
whole area contained within the half the radius of the 25 mag 
arcsec"^ isophote. The centre of the PPAK FOV is slightly dis- 
placed from the optical centre of the galaxy in order to include 
some interesting H 11 regions to the East of the centre. 

The PPAK fibre bundle consists of 331 science fibres of 
2.7" diameter each, covering a total hexagonal field of view 
(FOV) of 74" X 64" with a filling factor of ~ 60%. The sky 
is sampled by 36 additional fibres, distributed in 6 bundles of 



K. Viironen et al.: Spatially resolved properties of UGC 9837 



Table 1. Photometric properties of UGC9837 as observed by different surveys at different pass-bands. All the fluxes are given in 
Janskys. 



SDSS 


2MASS 


IRAS 


GALEX 


u 2.83E-03 
g 7.23E-03 
r 1.19E-02 
1.45E-02 
z 1.51E-02 


J 1.03E-02 
H 1.20E-02 
K 9.37E-03 


12 yum < 5.82E-02 
25 yum < 6.05E-02 
60yum 2.83E-01 
100 yum 8.97E-01 


FUV 2.73E-03 
NUV 3.29E-03 



6 fibres each, located following a circular distribution at about 
90" of the centre and at the edges of the central hexagon. The 
sky-fibres are distributed among the science ones in the pseudo- 
slit in order to have a good characterisation of the sky. In addition 
there are 15 fibres used for cahbration purposes. 

The V300 grating was used, covering the wavelength range 
~ 3700 - 7100A , resulting in an effective spectral resolution of 
~ 8 A (FWHM). Three dithered 600 s exposures of UGC9837 
were taken following a pattern of ARA= +1.15" and ADec= 
+0.78" in order to have a complete coverage of the science field 
of view and to increase the spatial resolution. An additional 120 
sec exposures of the spectrophotometric standard Hz 4 was ob- 
tained in order to perform the flux calibration. Also bias frames, 
arc lamp exposures and sky exposures were taken for data re- 
duction purposes. 



each particular pixel within the new grid the flux is assumed to 
be a weighted average of the n-adjacent ones: 



G(x,y) = ^WiF{xi,yi) 



Data reduction was carried out using R3D (ISanchezl2006l) . m 
combination with E3D (ISanchezll2061 . The reduction consists 
of the standard steps for fibre-fed integral-field spectroscopy, in- 
cluding the bias subtraction, cosmic ray rejection, and extrac- 
tion o f the individual spectra by a Gaussian-suppression tech- 
mque (ISanchezll2006h in order to minimise the effect of cross- 
talk between adjacent fibres. The ext racted spectra w ere stored 
in a row-stacked-spectrum (RSS) file ( ISanchezl2004l) and wave- 
length calibrated. Differences in the fibre-to-fibre transmission 
throughput were then corrected by comparing the wavelength- 
calibrated RSS science frames with the corresponding ones de- 
rived from the twilight sky exposures. Finally, the sky was sub- 
tracted and the spectra were flux calibrated. 

After reducing the three individual pointings, the corre- 
sponding RSS frames were combined into a sing le frame using 
standard procedures included in R 3D (see also ISanchez et alj 
|201 It iMirmol-Oueralto et alJ I20T1I) . First, each position table 
was combined in a single one by co-adding to each pointing the 
differential offsets recorded by the telescope. PMAS is equipped 
with a precise guiding system that allows to perform relative 
offsets with an accuracy of 0.1". Then, a 10" radius aperture 
spectrum is extracted from the expected centre of the galaxy in 
each of the individual pointings, taking into account the offsets. 
These three spectra are compared, and if a difference larger than 
a 5% is found, a rescaling is applied to match the intensity level 
of the first pointing. With this procedure possible variations in 
the transparency between the different exposures is compensated 
for. 

The resulting data were then spatially re-sampled in a regu- 
lar grid data-cube, by interpolating the information in each indi- 
vidual spatial pixel at each individual wavelength. Thi s is done 
by us ing the interpolation routines included in E3D jSanchezI 
I2OO4I) . by adopting a Natural-Neighbour non-linear interpolation 
scheme and a final pixel scale of T'/pixel for the resulting data- 
cube. The Natural-Neighbour routine adopted here is based on 
the Voronoi tessellation of a discrete set of spatial points. For 



1=0 

where the n-neighbours are selected as the adjacent ones if the 
new pixel was included in the previous grid and a Voronoi tes- 
sellation was performed. The weights correspond to the fraction 
of area that the new tessellation will remove from the old one, 
if the new pixel is included. These weights are renormalised in 
order to conserve the integrated flux, by taking into account the 
number of adjacent pixels considered in each interpolation. It 
is important to note here that this interpolation does not require 
to fill gaps between the fibres, which would be the case if only 
o ne single po inting procedure was adopted. As akeady shown 
in lSanchez et al. (2007), where this procedure was presented for 
the first time, the dithering guarantees a complete coverage of 
the FOV, and the interpolation routine resamples the data, with- 
out creating new values in areas not covered by the IFU. 

Finally, the flux calibration of the resulting datacube was ad- 
justed using the SDSS photometry. This was done because even 
though the relative flux calibration of our data should be reli- 
able, there are many uncertainties that can affect the goodness of 
the absolute spectrophotometric calibration. The photon-noise 
from the low surface brightness regions of the galaxy and the 
sky-background noise affect its quality. Also the sky subtrac- 
tion entails inaccuracies. In addition, the night of the observa- 
tions of UGC9837 was not photometric. Recalibrating our data 
based on the SDSS photometry available for UGC9837 is possi- 
ble, because our IFS data allows us to carry our spectrophotom- 
etry over a wide enough aperture so that the comparison with a 
corresponding imaging photometry can be carried out without 
worrying about the PSF effects. 

Of the five SDSS filters our spectrum covers the pass-bands 
of two, namely the g' and r' filters. We measured the counts in 
these SDSS images inside a 30" diameter aperture. Such aper- 
ture is large enough so as to avoid problems in the compared 
photometry regarding the exact position of the aperture centre 
and the exact radius of the aperture. These counts were con- 
verted to flux following the counts-to-magnitude prescription 
in SDSS documentatioiQ. The corresponding spectrophotome- 
try was then derived from our data, summing the flux of individ- 
ual fibres inside a 30" diameter circle and convolving this spec- 
trum with the SDSS g' and r' filter pass-bands. Using these two 
data pairs, an average scaling factor between the SDSS photom- 
etry and our spectrophotometry was calculated. This is shown 
in Fig|2] The error in SDSS photometry is estimated to be 0.15 
mag, including both the photon-noise and the background noise 
errors. The estimation is based on a statistical analysis in several 
different apertures of similar size on the back-ground subtracted 
8" X 8" post-stamp images. Our spectrum was then re-calibrated 



http://www.sdss.Org/dr7/algorithms/fluxcal.html#counts2mag 
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Fig. 2. The flux re-calibration based on SDSS photometry. Our 
spectrum of the 30" diameter central region of UGC9837, scaled 
by the average off'set between it and the SDSS g' and r' band 
photometry, is shown as a black line. The SDSS u',g', r', and 
band photometry are shown as blue dots and our spectrophotom- 
etry before re-scaling at g' and r' bands as green squares. 



using the resulting average scaling factor (0.94), and the new 
calibrated spectrum was used in the following data analysis. 



4. Data analysis 

The first step to extract any physical information from the data 
set is to identify the emission lines and to decouple their emis- 
sion from the underlying stellar population. After this, the emis- 
sion line spectrum is analysed using standard equations for 
ionised gas. 



4. 1 . Deriving tine emission iine fiuxes 

Population synthesis was used to model and subtract the stel- 
lar continuum. This technique results in emission-line measure- 
ments corrected (to a first-order) for stellar absorption. Many 
different techniques are now available to perform such a decou- 
pling. The one a pplied here is the one used in the analysis of the 
PINGS sample feosales-Ortega et alj|2010l) and the details of 
this process are described in Sanchez et al.l (1201 ih . Shortly, the 
scheme followed to decouple the stellar population and the emis- 
sion lines is the following: i) A set of emission lines is identified 
from the integrated spectrum, ii) For each spectrum in the data 
set, the underlying stellar population is fitted by a linear combi- 
nation of a grid of Single Stellar Populations (SSP), after correct- 
ing for the appropriate systemic velocity and velocity dispersion 
(including the instrumental dispersion which dominates the total 
observed dispersion when the grating V300 is used), and tak- 
ing into account the eff'ects of dust att enuation. The models were 
created using the MILES templates ( Vazdekis et alJl20T^ with 
three different ages (17.8, 1.0 and 0.09 Gyr) and two metallici- 
ties (Z = 0.004 and 0.03). Prior to the linear fitting, the nebular 
and sky emission lines are masked, iii) The fitted stellar popu- 
lation is subtracted from the original spectrum to get a residual 
pure emission-line spectrum, iv) Finally, the intensities for each 
detected emission line are derived. We note that the adopted SSP 
template library has an impact on the derived properties of the 
SSP, and actually is not suitable for reliably deriving the SSP 
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Fig. 3. Error in the derived Ha/H/? ratio as a function of the S/N 
ratio of the Hy6 line for a young (1 Gyr, Top) and old (17.78 Gyr, 
Bottom) SSP with three diff'erent dust values. Ay - 0.4 (red), 0.5 
(green), and 0.6 (blue). 



properties, but has little efi'ect on the der ived properties of the 
ionised gas (iMarmol-Oueralto et al.ll20TTl) . 

Next, the individual emission-line fluxes were measured 
in each spectrum by considering three spectral window re- 
gions: 1) [Oiii]3727A, 2) from Hy to [Oiii]5007A, and 3) from 
[Nii]6548A to [S n]6731A. The windows were selected so that 
they cover all the nebular emission lines detected in the inte- 
grated spectrum. A simultaneous multi-component fitting was 
performed using a single Gaussian function (for each emission 
line contained within each window) plus a low order polyno- 
mial (to describe the loca l continuum and to s implify the fitting 
procedure) using FIT3D (ISanchez et al.ll2066h . The central red- 
shifted wavelengths of the emission lines were fixed and since 
the FWHM is dominated by the instrumental resolution, the 
widths of all the lines were fixed as well. This procedure de- 
creases the number of free parameters and increases the accu- 
racy of the deblending process (when required). Line intensity 
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Table 2. Measured and dereddened emission line fluxes (Flux, Flux_dr, respectively) of the SDSS spectrum , and central and total 
integrated our spectra (SDSS, cent, tot, respectively) of UGC9837. All fluxes are first corrected for the underlying stellar absorption 
and then normalised to F(Hj6) = 1.0 and their errors are given in parenthesis. 



Line ID 


r [A] 


Flux (SDSS) 


Flux_dr (SDSS) 


Flux (cent) 


Flux.dr (cent) 


Flux (tot) 


Flux_dr (tot) 


[On] 


3726.8 






5.91 (0.09) 


13.67 (0.88) 


3.42 (0.03) 


4.83 (0.10) 


Hy 


4340.5 










0.43 (0.02) 


0.51 (0.02) 


US 


4862.0 


1.00 (0.06) 


1.00 (0.06) 


1.00 (0.06) 


1.00(0.06) 


1.00 (0.02) 


1.00(0.02) 


[Om] 


4959.5 


0.23 (0.02) 


0.22 (0.02) 


0.17 (0.02) 


0.15 (0.01) 


0.45 (0.01) 


0.44 (0.01) 


[Om] 


5007.8 


0.68 (0.05) 


0.65 (0.05) 


0.50 (0.04) 


0.45 (0.04) 


1.37 (0.02) 


1.31 (0.02) 


[Nil] 


6548.1 


0.55 (0.02) 


0.42 (0.03) 


0.65 (0.02) 


0.30 (0.02) 


0.24 (0.01) 


0.17(0.01) 


Ha 


6562.8 


3.71 (0.08) 


2.86 (0.19) 


6.21 (0.09) 


2.86 (0.17) 


3.93 (0.02) 


2.86 (0.05) 


[Nn] 


6583.6 


1.64 (0.06) 


1.26 (0.09) 


1.96 (0.07) 


0.90 (0.06) 


0.72 (0.02) 


0.52 (0.02) 


[Sii] 


6717.0 


1.19(0.07) 


0.90 (0.08) 


2.16(0.07) 


0.95 (0.07) 


0.91 (0.02) 


0.65 (0.02) 


[Sii] 


6731.0 


0.91 (0.05) 


0.69 (0.06) 


1.70 (0.07) 


0.74 (0.06) 


0.67 (0.02) 


0.48 (0.02) 



" Rest-frame 



Table 3. Physical properties of UGC9837 derived from the SDSS spectrum, and central and total integrated our spectra (SDSS, 
cent, tot, respectively). 



Property 


SDSS 


cent 


tot 


log([Nn]/HQ') 


-0.35±0.02 


-0.50±0.02 


-0.74+0.01 


log([Om]/H;e) 


-0.39±0.03 


-0.70±0.04 


0.312+0.006 


log(U) 




-4.21+0.04 


-3.473+0.009 




0.38±0.09 


1.13+0.08 


0.46+0.03 


K' [mag] 


0.82±0.20 


2.43+0.18 


1.00+0.06 


12 + log(0/H) (03N2, Pe 04) 


8.68±0.27 


8.68+0.27 


8.46+0.26 


12 + log(0/H) (S2N2, Vi 07) 


8.42±0.21 


8.25+0.19 


8.13+0.17 


12 + log(0/H) (Pi 10) 




8.36+0.16 


8.40+0.11 


SFR[Mg yr-'] 


0.06 


0.13 


2.48 




40 20 -20 



A RA (arcsec) 

Fig. 4. 2-D map of the signal-to-noise ratio of the HyS line. 
Contours correspond to the V-band emission derived from the 
datacube at the flux levels: 0.025, 0.035, 0.053, 0.077, 0.105, 
0.137, 0.172, 0.210, 0.251 and 0.295 x 10"^ erg s"' cm^^. 



fluxes were then measured by integrating the observed intensity 
of each line. 

For the error analysis, an average continuum RMS of the 
original spectrum was calculated and added by bootstrapping 
method to the SSP-Hgas model spectrum provided by the fitting 
program. This was repeated 100 times and each time the fitting 
was done on the resulting spectrum. The sigma of the median 
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Fig. 5. The N203 diagnostic diagram for our central spectrum 
(cent), the SDSS spectrum (sdss), the total integrated our spec- 
trum (int), and the individual fibre-spectra of signal-to-noise ra- 
tio > 5cr and Ha line intensity > 3 x 10"'^ erg s"' cm"^ arcsec"^ 
(blue dots). The solid and dashed lines separate the zones of star- 
burst gala xies (below th e lines ) an d AGNs (above the line s) as 
defined bv lKewlev etal] (1200 ll) and lKaufl^mann et al.l (l2003h . re- 
spectively. 



absolute deviation of these 100 simulations was then assigned as 
the eiTor in the measured emission/absorption line fluxes. 

To test if, and how much, the decoupling procedure used 
affects the derived emission line fluxes, the following experi- 
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Fig. 6. The integrated spectrum of UGC9837. The total inte- 
grated spectrum is shown in red, the best fitted combination of 
SSP models in black and the spectrum of the gas only in blue. 



ment was carried out. Different known emission line spectra (32 
spectra derived from different parts of M74) were summed to a 
m odel spectra of both y oung (1 Gyr) and old (17.78 Gyr) SSP 
of IVazdekis etaP (l20T0). assuming three values of dust, Av - 
0.4, 0.5, and 0.6. The resulting spectra was then decoupled us- 
ing the decoupling procedure introduced above, and the emission 
line fluxes measured. This was repeated 50 times adding random 
noise to the spectrum to get an idea about the involved errors. 
The derived emission line fluxes were then compared with the 
corresponding original fluxes. Fig. |3] shows the resulting differ- 
ence (original - measured) in the measured Ha/Hp flux ratio. We 
see that the errors in the case of the young SSP are somewhat 
bigger as compared to the old SSP, but in both cases the errors 
are reasonable, < 20% for S/N> 5. In the case of the young SSP 
one is also biased to under-estimate the H/3 emission flux due to 
underlying strong absorption. This bias can be expected to affect 
the extinction map derived in Sec. 15.21 from Ha/Hp ratio. The 
Ha line map should not be strongly affected due to the strength 
of the Ha emission, while the rest of the maps are not affected 
by the bias as line ratios are used in their derivation. We show 
in Fig.|4]the 2-D map of the signal to noise ratio of the H/S line, 
so that the reader can have an idea of the expected spatial dis- 
tribution of the errors associated to the decoupling procedure of 
ionised gas and stellar populations. In the center of the galaxy 
and along the arms the S/N ratio is generally > 5. 



4.2. Methods 

The emission line fluxes derived above reflect the physical prop- 
erties of the galaxy. Its dust extinction can be derived from the 
Ha/H/3 line rat io. Throughou t the pa per, we have adopted the ex- 
tinction law of lCardeUi et alj(ll989l) with Ry -3.l.We assumed 
the Balmer decrement for case B recombination of ionised gas of 
temperature 10^ K and density of 100 cm"^ (Osterbrock 1989), 
i.e., I(lia)/l(iy3) ~ 2.86, as the unextincted expected value. 

The nature of the ionisation is explored based on the 
Iog(rNn1/Ha) vsJog([0 iii]/Hy6) (N203) diagnostic diagram o f 
iBaldwin et all (Il98ll) (see also IVeilleux & OsterbrocM 119871) . 
This diagnostic is insensitive to dust extinction, as the needed 
emission lines are very nearby in wavelength. This diagram al- 
lows to separate the ionisation dominated by star-formation from 



the one dominated by an Active Galactic Nuc leus (AGN, see 
e.g.. lKewlev et aI.ll200lHkaufi^mann et ani2003l) . 

The degree of the ionisation is traced by the ionisa- 
tion parameter U. The lines available in our spectrum al- 
low us to derive the ionisation parameter from the ra- 
tio rOii1/fOm1^ A3727/(^4959 + ^5007) using the formula 
(iDfaz et al.ll20()0l) : 

logt/ = -0.80 Iog([On]/[0 III]) - 3.02. 

The electron density is derived from the 
[S n]6717A/6731A doublet ratio assuming a typical elec- 
tron temperature for star forming regions, = 10000 K (e.g. 
Osterbrock 1989) . 

For oxygen abundance determination we use two empirical 
calibrations: O3N2=log{([Om]^5007/Hy3)/([Nii]^6583/HQ')} in- 
dicator of Pettini & Pagel (2004) and the new calibration rela- 
tions of IPilvugin et al.l (12010 ) based on strong emission lines 
of oxygen, nitrogen and sulphur. The latter calibration relations 
are given separately to three classes of Hn regions: cool, warm 
and hot, and the criteria how to classify the H 11 regions is given 
in|Pilvugin et al. (2010). Even though we are not studying in- 
dividual Hi: regions here, we use the same criteria to decide 
which equation to use for deriving the oxygen abundance from 
each spectrum. Finally, we compare the results from the em- 
pirical calibrations to the ones derived from the model based 
log(Ha/[Sii]Xl6717+6731) vs. log(HQ'/[Nii]^6583), S2N2, in- 
dicator of IViironen et al.l (i2007). Ins tead of more commonl y 
used methods, such as the R23 index (iKewIev & Dopitall2002h . 
we choose to use the 03N2 and S2N2 indicators because the 
SDSS spectrum with which we are comparing our data does not 
cover the [On]/l3727 line. 

Finally, the star formation rate (SFR) of the galaxy is stud- 
ied based on its Ha line intensity. The extinction corrected line 
fluxes are transformed to absolute luminosities and the corre- 
sponding SFRs are then derived adopting the classic al relation 
between the SFR and the luminosity of the Ha line (iKennicuttI 
II998h . 

5. Results 

In this section the results of the data analysis introduced above 
carried out for integrated and spatially resolved spectra of 
UGC9837 are presented. 



5.1. Integrated spectra 

To create an integrated spectrum, all the individual fibre-spectra 
were convolved with the Johnson V filter pass-band and only 
those spectra with fluxes greater than 10""^ erg cm"^ s"' were 
summed. This was done in order to not to take into account the 
data of the fibres that do not contain much signal or do not con- 
tain signal at all as the fibres are sampling the regions where the 
intrinsic flux from the galaxy is low or null (e.g. the borders of 
the galaxy, intra-arm regions). A similar approach was adopted 
in t he analysis of t he inte grated spectrum of NGC628, presented 
by Sanchez et all (l20Ill within the PINGS survey. We derive 
two integrated spectra, one covering only the central parts of the 
galaxy and another covering the whole data cube and compare 
the physical properties derived from these spectra. For compari- 
son, we carry out similar analysis for the SDSS spectrum, which 
covers only a small central region of the galaxy. 

To derive the central spectrum, the entire pixels inside a 6" 
diameter are summed. The diameter of the SDSS spectrum aper- 
ture is 3". Since we are only interested in seeing quantitatively 
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the differences in the physical properties derived in the centre 
as compared to those derived when (almost) the whole galaxy 
is covered, we did not try to exactly match the SDSS spectrum 
aperture because it would require interpolation (and the corre- 
sponding errors) as our pixel size is 1". In addition, we do not 
know the exact position of the SDSS aperture. 

The total integrated spectrum of UGC9837 was created sum- 
ming the spectra of all the fibres (fulfilling the V-band flux crite- 
rion above). This spectrum, with the ionised and stellar compo- 
nents separated, is shown in Fig.|6] 

Once the integrated spectra are created, the stellar and 
ionised components are separated and emission lines measured 
and analysed as described in Sec. |4] Table |2] lists the emission 
lines and their fluxes with the corresponding errors in the SDSS 
spectrum, in our central spectrum, and in the total integrated 
spectrum. Table [3] lists the physical properties with their cor- 
responding errors derived for UGC9837 from the three above 
mentioned spectra. The errors in the oxygen abundance include, 
in addition to the the total errors propagated from the errors in 
our data, the errors in the models themselves: 0.25 dex for the 
03N2 method (Pettini & Pagel 2004), 0.075 dex for the new cal- 
ibrations of |Pi lvugin et al. (2010), and 0.15 dex as indicated by 
IViironen et al.l (i2007i) for the S2N2 diagnostics. In addition, the 
N203 diagram is shown in Fig.|5] The derived densities for the 
total integrated spectrum, central our spectrum, and the SDSS 
specti-um, respectively, are N^= 63+52, 144+125 and 107+173 
cm^-^ . The errors are large because these low densities are at the 
lower limit where the [S ii] ratio is hardly sensitive to the density 
anymore. 

A clear difference is encountered when comparing the lo- 
cation of the integrated total spectrum with the central spectra 
in the N203 diagram, the former pointing to higher excitation 
than the latter This is not surprising considering that the total 
integrated spectrum includes the high excitation star forming re- 
gions of the spiral arms. Despite this difference, all the three 
spectra are located in the zone of star-formation in this diagram. 
Also the ionisation parameter U indicates higher ionisation for 
the total spectrum as compared to the central one. Because the 
SDSS spectrum does not cover the [Oii]/l3727 line, the degree 
of ionisation could not be derived from the SDSS spectrum. The 
extinction coefficient. Ay, derived from the total integrated spec- 
trum is lower than the corresponding central value. The value 
derived from our central spectrum is about three times the value 
derived from the SDSS spectrum. In principle this could hap- 
pen if in the SSP/ionised gas decoupling of our lower resolu- 
tion data, the H/? absorption was underestimated. However, this 
would lead to a systematically higher values of all emission line 
fluxes (normalised to Hp) in our spectrum as compared to the 
SDSS spectrum. This is not the case, as seen in Table |2] Thus, 
we explain the difference in the two extinction values by short 
spatial scale variations of extinction, as the aperture sizes of the 
two spectra differ and we do not know the exact position of the 
SDSS spectrum as respect to our spectrum. In Sec. 15.21 we will 
show that there is significant spatial variations of extinction just 
near the centre of the galaxy. 

There are clear shifts in the oxygen abundances derived us- 
ing the different methods, the values derived from the S2N2 in- 
dicator being the smallest and the ones derived from the 03N2 
indicator the highest. These biases probably derive from the fact 
that in the integrated spectrum, the fluxes of different lines can 
come from physically different zones of the galaxy, with differ- 
ent ionising clusters and metallicities. Also, the S2N2 method, 
the one giving the lowest abundances, is not yet calibrated in a 
consistent manner with the other bright-line methods. Finally, 



the presen ce of diffuse ionised gas (DIG) can modify the final 
line ratios. 'Moustakas & Kennicut3 (l2006bh showed that a strong 
DIG would increase both the N2 and S2 ratios, but less so the 
higher ionisation 03 ratio. As a result, both S2N2 and 03N2 in- 
dicators would slightly overestimate the metallicity. However, as 
the integrated spectrum is dominated by bright H ii complexes, 
emitting strongly in 03, the effect of DIG in the abundance de- 
rived using 03N2 indicator is not necessarily significant (see, 
for example Garcia-Benito et al. 2010). In the new calibration 
relations of iPilvugin et al.i (i2010) the metallicity is calculated as 

12 + log(0/H) = ao + fliP+ fl2 log 03 

+ fla log(N2/02) + a4log(S2/02), 

where 

P= — — — and 02 = [On]^^3727 + 3729. 
03 + 02 

The effect of DIG on the relation is not discussed, but consider- 
ing that its emission causes S2 and N2 to increase while 02, an d 
especially 03, are less affected (iMoustakas & Kennicuttl2006bl) . 
probably the presence of DIG leads to a sUght overestimation of 
metallicity when this relation is used. 

When comparing the total and central spectra, the oxy- 
gen abundances of the form er derived from th e 03N2 
jPettini & Page]|2004 and S2N2 (IViironen et al.112007 ') diagnos- 
tics lead to a lower oxygen abundance as compared to the corre- 
sponding values derived from the central spectra. This could be 
expected, as negative metallicity gradients are generally present 
in spiral galaxies. H owever, the abundances derived using the 
iPilvugin et al.l (1201 Oh calibration relations do not show such a 
trend, giving roughly the same result for the two spectra. In all 
cases the derived abundance is lower than the solar oxygen abun- 
dance (12+log(0/H)=8.7, Scott et al 2002). 

We derive an integrated SFR of ~ 2.4 Mo/year It is of 
the order of the values derived for spiral galaxies and consis- 
tent with t he picture of an unperturbed galaxy evolution (eg., 
'Usui et al. 1998). As expected, being an additive property, the 
star-formation in the more reduced aperture of the central spec- 
tra is lower. When normalised by the aperture of each spectra, 
the SFR derived from the SDSS spectrum is actually the high- 
est and then decreases when the aperture increases. This can be 
understood as in the larger aperture spectra the spectra of high 
star forming spiral arms is mixed with the spectra of intra-arm 
regions with low level diffuse Ha emission. 

5.2. Two-dimensional distribution of tlie gas properties 

The real power of 2D spectroscopy resides in its capacity to map 
the physical properties of objects along their spatial extensions. 
We create the emission line maps for UGC 9837 from the final 
reduced datacube carrying out the fitting procedure described 
in Sec. 14. II for each individual spectra. For each map, the pix- 
els with a flux < 0.2 x 10"'^ erg cm"^ s"' at the corresponding 
wavelength are masked, in order to retain only the high signal- 
to-noise data (> 5cr). When creating maps of emission line ra- 
tios, only those locations are considered where all the involved 
lines are above this selected minimum flux. Practically a map 
of any detected emission line, or a combination of them, could 
be derived from the individual fibre spectra. Here we present a 
number of maps corresponding to the most important emission 
lines and derived quantities of interest. 
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Fig. 7. Top Left: Ha line intensity map in units of 10 erg s ' cm"^ arcsec"^. Top Right: Dust extinction map of the ionised gas. 
Middle left: log([Nn]/HQ;) and Middle right: log(rOm]/H^) d i agnos tic line-ratio maps. Bottom Left: Dist ribution of the oxyge n 
abundance according to the 03N2 indicator of iPettini & Pagell (|2004 and Bottom Right: S2N2 relation of ' Viironen et all (l2007h . 
The contours represent the V-band image derived from the data-cube and show the flux levels 0.025, 0.035, 0.053, 0.077, 0.105, 
0.137, 0.172, 0.210, 0.251, and 0.295 xlO""' erg s"' cm^l 
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As the dominant gas ionisation mechanism giving rise to the 
emission Hnes in this galaxy is due to hot (OB) stars, the Ha line 
intensity map, shown in the top-left panel of Fig. |7] traces the 
star forming regions of the galaxy. The map shows a clear spiral 
structure with H ii complexes of different sizes and morphology 
along the spiral arms. The brightest sources are located in the 
Northern arms. Also a clear presence of diffuse emission along 
the spiral arms and in the intra-arms region is visible. 

The dust extinction was derived from HafH/3 line ratio, as in 
the previous section, and is shown in Top right of Fig. [T] in terms 
of Ay- The dust extinction is very variable, but as a general trend 
it is strongest in the bulge and then follows the spiral arms. We 
note that in the areas populated by young stars, i.e. expectedly in 
the spiral arms, the extinction may be overestimated by < 20% 
(see Figs. [3] and |4]) A strong short-scale variation of extinction 
in the centre of the galaxy is visible (see also Figure [D, as was 
already indicated in Sec. l5.ll bv the differences between the ex- 
tinction values derived from our and SDSS central spectra. 

The dust extinction map described above is used to deredden 
the rest of the line maps and only the pixels covered by the dust 
map, in addition to the data on the lines in question, are consid- 
ered. The middle panels of Fig. |7] shows t he diagnostic line-ratio 
maps log([Nn]/HQ'), and log([Om]/H/?) dVeilleux & Osterbrockl 
[T987h . The radial gradient is clear, the ionisation hardening to- 
wards the edge. This is indicated by the strengthening of the 
[O in]/Hj6 ratio towards the edge while [N ii]/Ha ratio becomes 
smaller No clear trend along the spiral arms is visible. As 
shown in Fig. |5] the values of both ratios at any Ha emitting 
location in the galaxy are consistent with ionisation pro duced 
by hot OB stars according to both iKewlev et alJ (l200il) and 
iKauffmann et~aD f2003). 

Next we studied the spatial distribution of the oxygen abun- 
dance as for the integrated spectra above using for brevity only 
the 03N2 indicator of Pettini & Pagel (2004) and S2N2 diag- 
nostics of Vi ironen et alj (I2007D . The resulting maps are shown 
in the bottom panels of Fig. |7] As the comparison of the cen- 
tral and integrated spectra above already hinted, in both maps 
the oxygen abundance is clearly strongest in the centre and de- 
creases towards the edge, showing a gradient of ~ -0.4 dex, but 
03N2 about 0.3 dex offset from S2N2. 



5.3. Simulations of observations at liiglier redsliifts 

In this section we study how the above maps would 

look like if the galaxy would be located at a higher 

redshift (z - 1-2) and o bserved as part of sur- 
veys such as SINS jpo rster Schrei ber et al 



MASSIV (EDinat etaDl2009t lOuevrel et al 



2006) or ongoing 
2009), both using 



VLT/SINF ONl. For this we have u sed a code adapted from the 
one used in lEikenberrv et all (l2006h . The simulated spatial scale 
is 0.125"/pix, and the data was depixelised accordingly, sum- 
ming only entire pixels. Considering the original and simulated 
pixel sizes there is no differences in the depixelisation for sim- 
ulated redshifts of 1.5 (typical for MASSIV) or 2.2 (typical for 
SINS). Next the flux was rescaled as 



^simul 



l+z. 



■orig 



1 Zsimul I 



where dorig and t/i,,,,,,/ are the original and simulated luminosity 
distances, respectively, and Zorig and Zsimui the original and sim- 
ulated redshifts. Finally, we assumed seeing conditions of 0.6" 
and convolved the data with two-dimensional Gaussian kernel 
of 0.6" FWHM. We accepted for the resulting datacube only the 
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Fig. 8. Closer look at the extinction variation (see Fig. QTop 
Right) in the centre of the galaxy. 



pixels where the convolution does not rely on the zero-padding 
used outside the image area. We have not added any noise to 
the simulated data as it depends on the observation strategy (e.g. 
integration time). However, in real observations at high redshift 
one would generally expect lower S/N ratios than in the lower-z 
observations. We neither considered the instrumental response 
of the simulated instrument at the wavelength where the line 
fluxes would be shifted at the simulated redshift nor we consid- 
ered if the galaxy is intrinsically bright enough to be detected by 
the instrument at this redshift. Thus, to summarise, our simula- 
tion shows what kind of data could be obtained on a galaxy like 
UGC9837 if it was detected by an IFU of 0.125" pixel scale in 
typical observing conditions of 0.6" seeing. At the level of this 
simulation the picture would change only slightly if instrument 
like Keck/OSIRIS with 0.1" pixel scale was simulated while if 
a smaller pixel scale (which generally leads to a smaller FOV) 
was assumed, the result would be improved. 

In Fig. |9] we show the simulated maps (at z-2.2) derived 
from the maps in Fig. |7] As there would be no differences in 
the depixelisation and as we have not touched the S/N ratio, nor 
considered the instrumental response, the corresponding maps 
at z= 1 .5 would be equal but of higher flux. Comparing the maps 
in Figs. |2]and|9] we see that in the high-redshift simulation the 
spatial information about the SFR (Ha line map) is lost except 
for a gradient of stronger Ha emission towards North, where the 
brightest H ii regions are located. The information about the dis- 
tribution of dust is lost completely. 

The combined effect of the latter two not only prevents us 
from deriving information about the SFR distribution along the 
galactic disc, but also lead to a bias in the level of SFR itself. As 
discussed in l5.2l the spatial distribution of the extinction by dust 
is very variable, being strongest in the centre and then following 
the spiral arms. The loss of this information leads to underesti- 
mation of absolute line fluxes, and SFR, in the dust rich parts of 
the galaxy, while opposite is true for the less dusty parts. 

As for the relative line fluxes, the bluer line intensities are 
underestimated as respect to the redder lines in the dusty regions, 
especially in the centre and in the spiral arms. For this reason, we 
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Fig. 9. The same maps (but in units of 0.5 x 10 ^'^ erg s ' cm ^ arcsec ^) as in Fig|7]but as they would be seen by VLT/SINFONI if 
the UGC9837 was placed at redshift z=2.2. 
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warn against using dust dependent indicators, such as emission 
line ratios of lines with large wavelength difference, or absolute 
line fluxes, if the distribution of dust extinction is not known. 

The gradients in log([Nii]/Ha) and log([Oiii]/Hj6) line ratio 
maps as well as in the 12 + log(0/H) can be still roughly appreci- 
ated in the simulated maps. However, the gradients appear shal- 
lower in the simulated maps as compared to the original ones. 
Shallower gradients can be expected as the spatial resolution of 
the simulated maps is worse so that zones of higher and lower 
ionisation and metallicity are mixed. This has its effect in the 
interpretation of the abundance gradients derived from observa- 
tions with limited spatial resolution. 

The inside-out growing scenario of galactic discs predicts 
the presence of negative metallicity gradients in spiral galax- 
ies. However, different chemical evolution models disagree on 
the time variation of these gradients. Some models predict that 
the gradients steepen with time (e.g. Chiappini et al. 2001) while 
other models predict the opposite (e.g. Fu et al...2009 ). Recently, 
lOuevrel et al.l (l201 ih published a study of metallicity gradients 
in a sample of galaxies at z ~ 1 .2 based on SINFONI data. They 
found both negative, zero and positive gradients, explaining the 
latter by interactions. If a gradient was detected, it was in all 
cases shallow as compared to typical gradients of the nearby 
galaxies. Based on our simulation above, we argue that the gra- 
dients they measure are rather lower limits of the real abundance 
gradients of these galaxies. 



6. Summary and conclusions 

We have caiTied out a detailed study of the ionised gas 
of UGC9837 using the 2D data obtained with PMAS/PPAK 
dMarmol-Oueralto et alj|20^1 1,) . We have studied the integrated 
and spatially resolved properties of the ionised gas, deriving the 
nature of ionisation, ionisation parameter, oxygen abundance, 
SFR, and dust extinction (Ay) of the galaxy. First we compared 
the ionised gas properties derived from a spectrum covering only 
the centre of the galaxy (derived from our data and a SDSS spec- 
trum) with the corresponding properties derived from a spec- 
trum covering (almost) the whole galaxy. Clear differences were 
encountered as the spectra are sampling different zones of the 
galaxy. We conclude that deriving global properties of the galaxy 
from a spectrum of only partial spatial coverage would lead to 
significant biases. We derived for the galaxy a total SFR of ~ 2.4 
Mo/year which is consistent with the picture of an u nperturbed 
galaxy evolution for this object (eg.. lUsui et alj|T998i) . The in- 
tegrated oxygen abundance derived using different strong line 
methods varied but all diagnostics lead to a subsolar abundance. 

Next, the spatially resolved properties of the ionised gas 
were studied. We found that the ionisation becomes harder with 
increasing distance from the centre of the galaxy, while negative 
gradients of both oxygen abundance and dust extinction where 
uncovered. The Ay distribution is consistent with tracing spiral 
arms, except for one (maybe two) arm(s), for which the dust 
distribution points to a leading nature of the arm(s). In general, 
our results are consistent with an inside-out growing scenario 
for disk galaxies, where the more evolved stellar populations 
and richer gas components with more dust, are located in the 
centre of the galaxies, and the less evolved stellar populations 
and metal poor gas, with less dust grains and stronger SFR and 
harder ionisation is located in the outer regions. 

We also carried out a simulation how the galaxy would be 
seen if it was placed at a higher redshift and observed as part of 
surveys such as MASSIV or SINS. In our simulation the spatial 



information about the SFR and dust distribution is lost while ra- 
dial gradients of ionisation strength and oxygen abundance can 
still be roughly distinguished. The gradients derived from the 
simulated data are shallower than the original (more realistic) 
gradients. These biases should be kept in mind when interpreting 
observations with limited spatial resolution. According to our 
Local Universe data, the dust distribution along the galactic disc 
can be very variable. If this information is lacking, dust inde- 
pendent indicators should be used for the study of the galaxy in 
question. The metallicity gradients derived from data with poor 
spatial resolution, in turn, are lower limits of the actual ones. 

Finally, we want to remark the importance of a statistical 
sampl e of 2D data on nearby galaxie s. The already existing data 
of the lMarmol-Oueralto et al. (|20I Ih sample, and ultimately the 
data on a large sample of Lo cal Unive rse galaxies being obtained 
by ongoing CALIFA (Sanch ez et ^1 12010: Sanchez et al. 201 1]) 
survey will allow analysis like the one presented here to be car- 
ried out in a statistical manner. This in turn will provide strong 
constraints on the theories of galaxy evolution and fix observa- 
tional properties of galaxies in the Local Universe, which may 
have a potential impact in the interpretation of observed proper- 
ties at higher redshifts. 
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